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A robust scheme for state preparation and state trapping for the spins of two electrons in a 
semiconductor double quantum dot is presented. The system is modeled by two spins coupled to 
two independent bosonic reservoirs. Decoherence effects due to this environment are minimized 
by application of optimized control fields which make the target state to the ground state of the 
isolated driven spin system. We show that stable spin entanglement with respect to pure dephasing 
is possible. Specifically, we demonstrate state trapping in a maximally entangled state (Bell state) 
in the presence of decoherence. 
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I. INTRODUCTION 

In recent years, there has been a great interest in quan- 
tum computing and quantum information processing. 
Since Shor demonstrated that his quantum algorithm is 
able to perform prime factorization more efficiently than 
any known classical algorithm, the attention paid to this 
type of computation has grown enormously! 1 ^ 

Within a network model, the fundamental building 
block of a quantum computer is the quantum bit or qubit. 
A two-level system (TLS) which can provide coherent 
state superposition is needed to physically realize a quan- 
tum bit. There are several suggestions in the literature 
for achieving the goal of finding a suitable TLS. Amongst 
such proposed systems are quantum dotsp' single pho- 
tonsPthe nuclear spinP superconducting devicesP^and 
trapped ionsP 

Solid state realizations of qubits and quantum gates 
are attractive because they may be quite easily scaled 
to larger arrays of qubits which is necessary for efficient 
quantum computing. A major advantage of the semi- 
conductor solid state implementation is that the existing 
semiconductor technology is well developed. Because of 
these advantages, we are going to present our approach 
in terms of semiconductor spin quantum dots, although 
the model and strategies developed here can be mapped 
onto most qubit realizations which have been discussed 
in the literature. 

Quantum computation is powerful because it takes ad- 
vantage of quantum interference (i.e., coherence) and 
entanglement. However, when a quantum system is in 
contact with an environment, the resulting interaction 
destroys coherence and changes the populations. Sev- 
eral methods for suppression or avoidance of (unwanted 
effe cts of) decoherence, such as quantum error correc- 
tion) 1 ^ encod ing of the qubit into a d ecohe rence free 
subspacej 10 * 11 ! and dynamical decoupling] 12 * 13 ! have been 
proposed. Another technique to cope with dissipation 
is to make the desired state to the ground state of the 
system by tuning of (possibly adiabatic) control fields. 14 

In this paper, we study the feasibility of Bell-state gen- 
eration for two electron spins in a semiconductor double 



quantum dot. The present work is structured as follows. 
In the first part, we outline the model which is used to 
describe the system of interest. The second part is de- 
voted to the derivation of the equations of motion and 
the calculation of the system-bath correlation functions. 
In the next part, the strategy for fighting decoherence 
is outlined and numerical results are given. We will end 
with summary and conclusions of the present work. 



II. THEORY 

The spin \ of an electron provides a natural two-level 
system. An electron in a semiconductor quantum dot 
as a realization of a qubit has been proposed by Loss 
and DiVincenzoP The spin directions up and down with 
respect to an external magnetic field represent the two 
basis states of the qubit. Few-electron semiconductor 
quantum dots can be realized, for example, by means of 
surface gates on top of a GaAs/AlGaAs heterostructure 
which holds a two dimensional electron gasP^ The num- 
ber of conducting electrons in each dot can be controlled 
and monitoredP21 

It has been shown that a system consisting of two 
neighboring quantum dots, each populated by one excess 
electron, can be described by a Heisenberg-like Hamil- 
tonian H s ' = J(t)Si ■ 5*2 if the dots are coupled via a 
tunable tunneling barrier plfil The straightforward way to 
implement such a two-qubit array is to use a heterostruc- 
ture as described above and to increase the number of 
gate electrodes PS3 in a recent paper it has been argued 
that for quantum dots (which are spatially separated so 
that tunneling between the dots can be neglected), con- 
trolled entanglement via electrostatic interaction is pos- 
sible.^ 

First, the Hamilton operators for the double quantum 
dot, the environment, and the interactions on which our 
investigation is based are defined. To include the influ- 
ence of the environment in the dynamics of the double 
dot, a Mar kovian quantum master equation approach is 
used below P^ 19 l 2 ° lThe system-bath correlation functions 
are calculated analytically. 
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A. Description of the system 

The model system consists of two conduction electrons, 
each sitting in its own quantum dot. Each is described 
by a Hamiltonian of the form 



H<P(t) 



-g* —S z B?{t) 
2m* 



9_ 
' 2 



(1) 

(2) 



g* denotes the gyromagnetic ratio which depends, as does 
m*, on the type of the semiconductor quantum dot. This 
Hamilton operator describes the interaction of the ith 
electron spin with an external magnetic field applied in 
the z direction. B^\t) is a control field which is used to 
adjust the Zeeman splitting associated with the electron 
spin in quantum dot i. For the most general case of 
magnetic fields, there are nonzero x and y components, 
e.g., to perform single qubit operations on the spin. 

We model the presence of an environment which is re- 
sponsible for dephasing by coupling the electron spin de- 
gree of freedom to a bath of bosons. Such models have 
been widely used in the literature to describe various 
forms of interactions at a microscopic level by selection of 
the spectral properties of the bosonic bath. In the present 
system, the main cause for dephasing arises from charge 
fluctuations in the vicinity of the q uantum dots , phonons, 
and interaction with nuclear spinsP ^ 2 * 2 ^ 24 * 25 ! The influ- 
ence of the phonons on the spin degree of freedom is me- 
diated through the spin-orbit interaction. At low tem- 
perature, acoustic electron-phonon interaction becomes 
the main scattering mechanism! 2 ^ Therefore, we can as- 
sume that only LA and TA phonons are involved in the 
scattering process. As all nuclear spins corresponding 
to III-V semiconductor materials, such as GaAs, have 
nonzero nuclear spins, the hyperfine interaction of the 
latter with the electron spin is another source of decoher- 
ence.^ Moreover, for exchange-interaction-based quan- 
tum gates, charge fluctuations produced by tuning the 
gate voltages lead to dephasing.^ 

We model the main effects from these interactions, i.e., 
dephasing, using two uncorrelated baths of harmonic os- 
cillators, 



H; 



1,2, 



(3) 



where 6^ and b't 1 are the bosonic creation and anni- 
hilation operators for modes with frequency . The 
inter action of the baths with the qubits is of spin-boson 



type, 1 



27 28 



(4) 
(5) 



(i) 

The g k are the effective coupling constants of the spin- 
boson interaction. For the case of strongly correlated 
reservoirs (collective decoherence), the decoherence free 
subspace includes a set of maximally entangled states 
between the two spins (see Sec. II B I. In Ref. [24l it has 
been shown that for the case of diagonal qubit -qubit cou- 
pling and collective decoherence a set of robust entangled 
states can be identified. So, here we examine the more 
challenging (and probably more realistic) case of two in- 
dependent reservoirs. In addition, we are going to im- 
plement the full Hcisenberg interaction between the two 
spins. 

The coupling to in Eq. Q leads to pure dephas- 
ing, i.e., destruction of coherence. Terms proportional 



to a x or <jy' would cause population relaxation which, 
within this model, corresponds to spin flips. Spin flips 
do occur in semiconductors. However, these processes 
in general take place on a much larger time scale than 
dephasing processes, i.e., the spin relaxation time T\ is 
much longer than the dephasing time T%, and so relax- 
ation is neglected here. In this sense we study the quan- 
tum dots on a time scale t <T\. 

The qubit-qubit interaction which is needed to pro- 
duce entangle ment and conditional operations is of 
Heisenberg typep^ 



(6) 



a« • a& = o-W ® 4 2 ) + 4 1 ) ® 4 2 ) + ® af\ 



The Hamiltonian of the externally controlled two-qubit 
system is 

H s (t) = H s (1 \t) + H s (2) (t) + H s (12) (t) (7) 
= -fcBW(t)o-W ® l _ hBf\t)t ® + J(t)a^ ■ a (2 \ 

The system-reservoir Hamiltonian to introduce the inter- 
action between the environment and the system reads 

TT _ rr(l) , rr(2) 
HSR - H SR + H SR 

= ft (a« ® i) rw + £ (l ® erf )r< 2 \ (8) 

and the reservoir Hamiltonian consisting of the two in- 
dependent baths is of the form 



H R = H£' + H 



= E^W+E^rW'- (9) 

k k 

The full Hamiltonian of the interacting system is 



H(t) = H s (t) + H R + H S R. 



(10) 



The representation of H(t) in the system space is of 
block-diagonal form and we separate it into an "outer" 
and an "inner" contribution, respectively, H a (t) and 
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H b (t) 



H(t) = H a {t) + H b {t) 



"JTiiW 








H 44 (t). 





H 22 (t) 

Zf 2 * 3 (t) 





0" 

H 23 (t) 

H 33 (t) 





(11) 



With |0) and |1), respectively, denoting spins down and 
up, we use the following convention for the four basis 
states \i) ® \ j) = \ij), for i,j = 0, 1: 

|1) = |11> = (1,0,0,0) T , |2) = |10) = (0,1,0,0) T , 
|3) = |01> = (0, 0, 1, 0) T , |4) = |00) = (0, 0, 0, 1) T , 

i,j'=l,2,3,4 

for any operator A defined on the Hilbert space of the 
double-qubit system. One then finds 

H n (t) = -j5i(t) - E 2 (t) + j(t) + nrw + nrw + h r , 

H 22 (t) = -E x {t) + E 2 {t) - J{t) + hT^ - KT<® + H R , 
H 33 (t) = Exit) - E 2 (t) - J(t) - HT^ + nr® + H R , 
H u (t) = E x {t) + E 2 (t) + J(t) - ftr (1) - hT {2) + H 
H 23 (t) = H 32 (t) = 2 J, 



(12) 



pW, T( 2 \ and Hr are operators defined over the boson 
Hilbert space. 

In the latter part of this paper, we will use the den- 
sity matrix formalism of quantum mechanics. Density 
operators corresponding to the double-qubit system will 
be denoted by ps, whereas density matrices p without 
subscript S refer to the composite system (or subspaces 
thereof) including the bosonic environment. 



Since 



B. Equations of motion 



[H a {t),H b {t')\ =0, 



(13) 



the propagator for the composite system may be written 
in the form 



U(t,0) = U a (t)U b (t), 



(14) 



with 



-i f Hj(t')dt' 

Uj{t)=Te 8 , for j = a,b , (15) 



and T representing the time-ordering operator. It is 
instructive to define two special density matrices corre- 
sponding to the composite system for t — 0, 



P (a) (0) = 



pii(Q) p 14 (0)' 













AM Mo). 



o o 



0' 



p 22 (0) p 23 (0) 
Phi®) P3 3 (0) o 




(16) 



(17) 



similar to H a (t) and H b (t). One can calculate their time 
evolution as 



P ^\t) = U b {t)p^{Q)ul{t), 



(18) 
(19) 



by using 



H b (t),p^(0) = H a (t),pV»(0) =0. (20) 



('■')/ 



Thus, if we start with a density matrix of the form of 



Eq. (16) or (17 1 the system will remain in the corre- 



sponding subspace and we are effectively dealing with 
two noninteracting two-level systems. By applying the 
basis transformation given in Table [i] we are able to map 
the Hamiltonians H a (t) and H b (t) onto 



H r a ed (t) = -(E 1 (t)+E 2 (t))a z + J(t)l 

+ a, (firw+arw) +h r , 

H r b ed (t) = 2J(t)a z - J(t)t 

+ a x (nrw - hr^ + E 2 (t) - E^t) 



(21) 
(22) 



for one and the 



It can be seen that for I^ 1 ' = , 
same bath for both spins, the bath terms in Eq. (22 1 can- 
cel. Therefore, the "interaction part" of the Hamiltonian 
jjred vanjghgg ano ; n0 entanglement with the environment 
occurs. So, for the case of collective decoherence, density 
matrices of the form 



Ps 



PS, 11 









PS,22 

Ps.23 






PS, 23 
PS, 33 









PS,44 



(23) 



(PS 



tr R {p}) 



are decoherence-free^l even if we are dealing with a full 
Heisenberg-like qubit-qubit coupling (see Eq. (|6|). To 
be precise, the density matrix elements p5.11 and 
should be called stationary. This can be seen if we set 



Pu(0) = P4i(0) = in Eq. (16). Because the propagator 
U a is diagonal, it now follows that 



Ps,u(t) 

PS,iA{t) 



PS,ll{Q), 

Psm(0). 
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Furthermore, control in the p^ subspace via the qubit- 
qubit coupling is not possible because there is only a 
trivial (i.e., proportional to 1) dependence on J(t) in 
Eq. pi). 



Four-level system 


Reduced TLS 


Subspace 


|^+> = -^(0, 1,1, Of - 


(FOf 


P W 


|VO = ^5 (0,1,-1, Of - 


(o,i) T 




|11) = (i,o,o,of 


(1,0) T 


p (a) 


|00) = (0,0,0,1) T 


(o,i) T 





TABLE I: Mapping of the basis for the full four-level system 
to one for the reduced TLSs. 



If the initial state of the double dot is any valid super- 
position of j(/ a )(0) and p( b \0), it is possible to formulate 
and solve the master equations for each of the two sub- 
systems independently. However, in the following, we are 
going to solve the full four level dynamics. 

In the interaction picture (operators denoted by a 
tilde), using 



-i f (H s (t')+H R )dt' 
U o (t,0) = Te « l ' , 

p(t) = Ut(t,0)p(t)U o (t,0), 

H SR (t) = U$(t,0)HsR(t)U o (t,0), 



(24) 
(25) 
(26) 

the master equation in the Born-Markov approximation 
for the present system-bath interaction is of the forrrP 

f t Ps(t) = (27) 



1 



dt'tiR { [H SR (t), [H SR (t'),p s (t) ® p R (0)]] }■ 

When evaluating the double commutator one encounters 
the correlation functions 

(f «(t)f (i) (t')) R = tr R {f W(t)f (l) (t')p R (0)} . (28) 

These can be calculated analytically if we use an Ohmic 
spectral densityp^ES 



oc 



duj {J (lu) ...}, 



r\Loe~ 



(29) 



lu c is a cutoff frequency and 77 is a parameter which de- 
scribes the coupling strength of the bosons to the qubit. 
The specific value for lo c depends on the physical nature 
of the dephasing mechanism. 

Up to second order in the system-reservoir interaction, 
we may use the equilibrium form for p R in the calculation 
of the correlation functions, 



Pr 



Pit 



m 



Pr 



1 — exp 



(2) 

>Pr ■ 



With (3 — -jt^j, , we get 

( ri "<" f<a < ! ')> a ^{*'( 1 + 

, fl + iui c (t-f)' 



1 - iuj c (t - t') 



fiio c l3 



(31) 



where ip' is the derivative of the digamma functionPThe 
calculation of Eq. (31 ) is outlined in Appendix |A| 

Depending on the physical nature of the interaction, 
one may use more sophisticated formulations of the spec- 
tral density^ In Ref. |3TJ the dissipative interaction be- 
tween the orbital degrees of freedom and an acoustic 
phonon bath was investigated. This interaction leads to 
an effective spin-phonon interaction which is related to a 
more complicated spectral density. However, using these 
results would give correlation functions which cannot be 
calculated analytically. 

In the following calculations we will choose T — 50 mK 
for the temperature. At this temperature, the spin-flip 
energy AEz (at B = 1 T) is approximately five times 
larger than ksT. So, we can neglect spin flips and con- 
centrate on dephasing processes. Dephasing is an inter- 
esting phenomenon because decoherence associated with 
this process can occur even when there is no real energy 
exchange with the environment. Also, qu antum vacuum 
fluctuations can cause dephasing! 32 ! 33 * 34 -! 



C. Optimization strategy 

Assume first that we want to trap the system in the 
Bell state where 

|^ ± )-^(|1)®|0)±|0)®|1))= ~(0,1,±1,0) T . (32) 

Based on the discussion above, the strategy is to make 
this state to the nondegenerate ground state of the iso- 
lated driven (i.e., decoupled from the bath) spin sys- 
temPFrom Eq. ([22) and Table [fl we see that for J(i) < 



0, the state corresponding to indeed becomes the 

ground state of the nondissipative reduced TLS b. Thus, 
if we set J(t) large and negative, we would infer that this 
state is rather robust with respect to dephasing because 
the levelsplitting between and IVO is 4| J(t)\. (We 
set Ei(t) — E 2 {t) for simplicity.) 

This statement is readily verified numerically. Suppose 
that we start with the 2x2 density matrix for TLS b (see 
Sec.lilB), 



P$\0) 



1 




IV>+}(^ + | 



Without the qubit-qubit interaction (i.e., J 
evolves into the mixture 



(33) 
0), it 



(30) 



pf(oc) 



I 
5 



1 



1 
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Now, let us set J(t) = J = const and examine the 
dependence of the time-averaged purity, 



P(tf) 



1 



tr : 



[pf\tf)dt, (35) 



and the time-averaged fidelity. 



i 



F(t f ) 




P f(0) 



(36) 

on the qubit-qubit coupling J and the inverse temper- 
ature (3 SC = {huj sc = 1 meV) with Pg\o) given by 



Eq. (33 1. (In the following, physical quantities with a 
subscript sc are scaled with respect to uj sc .) The result 
can be seen in Fig. [T] keeping in mind that for bipartite 
systems, 



(37) 

15 K), good re- 
— 1 meV (qubit- 



= 1, 
1 

2' 

For "low" temperature (up to T 
suits can be obtained by setting J 
qubit couplings of several 100 /ieV have been reported in 
Ref. [15j). However, it should be noted that for T « 15 K 
and B z = IT, spin-flip processes constitute the main 
mechanism of decoherence and our model accounting for 
pure dephasing only is no longer valid. 

The conclusion of Sees. |IIB| and |IIC| may be sum- 
marized as follows. Given a Hamiltonian of the form 



of Eq. ( 12 1 , suppression of pure dephasing of the maxi- 
mally entangled states is possible by tuning of the 
qubit-qubit interaction. Explicit results of this strategy 
for state trapping and steering will be given in the next 
section. 

On the other hand, dephasing of states, 

I**) = -^(|1)®|1)±|0)®|0)) = ±= (1,0,0,±1) T , (38) 

cannot be reduced by means of changing J. To do so, the 
qubit-qubit Hamiltonian has to display anisotropy, e.g., 

H s ^ = J{t)a^®af\ (39) 

With the isotropic Hamiltonian [Eq. (fl2|)], it is also 
impossible to prepare dephasing-insensitive entangled 
states of the form 

1 ,„> . .... . 1 



. . ,,J)±i|0>®|l))= (0,l,±z,0) J 



1 



For this task, a coupling, 



(40) 



(41) 



would be needed. Hence, to trap an arbitrary Bell state 
by this strategy, a detailed control of the spin-spin cou- 
pling is required. Realization feasibility depends on the 
physical nature of the double qubit. 



In (Ac) 




meV 



(b) 



In (ftc) 




FIG. 1: (a) Time-averaged purity P and (b) fidelity F vs 
qubit-qubit coupling J and the logarithm of inverse temper- 
ature, n = 0.1, tf — -j- and lu c = 20 meV. 



III. NUMERICAL RESULTS 



We now investigate feasibility of the strategy discussed 
above at the example of state trapping and steering into 
the Bell state by solving Eq. (27 1 numerically. In 



order to quantify the degree of success, we use several 
well established quantities. A convenient measure for 
the degree of entanglement of bipartite systems is the 
concurrence.^ The concurrence C equals 1 if the qubits 
are maximally entangled. It becomes if the system 
factorizes or, for the case of a mixed density matrix, if the 
system can be represented by a mixture of factorizable 
pure states. Another important measure is the purity P, 
which gives information whether the system is in a pure 
state or in a mixtureP^ A purity of 1 indicates a pure 
state, while | < P < 1 occurs if the system is in a mixed 
state. 
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A. Trapping 

The aim is to steer the double-spin system into the 
entangled state and to trap it there starting from 

the initial state, 



(a) 



1^X^/1 = 





10 







IV> + XV> + I = 







(42) 

As outlined in the last section, we use variation of the 
qubit-qubit coupling J(t) to steer the system into the en- 
tangled state and to trap it by making it insensitive 
to pure dephasing arising from the interaction [Eq. Q]. 
We start with a control field J(0) = 0. The selected time 
dependence of the control field is shown in Fig. [27 f ) . Fig- 



V 


htv ac 


T tf = tf, sc /tjj sc 


huic 


bP 


BP 


0.08 


lmeV 


50 mK 13.04 ps 


20 meV 


1 T 


1 T 



TABLE II: Physical quantities used to calculate the data 
shown in Figs. [^a)-^!). 



ures. [2]ja)-[2jd) show that steering to and trapping in the 
entangled state |?/> + ) using this technique works very well. 
Convergence of the off-diagonal density matrix elements 
is generally better than for the diagonal elements. The 
concurrence shown in Fig. [2^e) demonstrates that prac- 
tically perfect entanglement (C m 1) can be achieved. 
Furthermore, the achieved entanglement is stable with 
respect to pure dephasing. 



B. Optimized driving and trapping 

In this section, we utilize a conjugate gradient (CG) 
algorithm to optimize the control field plotted in 
Fig. j3^e)F^ This can be done by formulation of a cost 
functional 



T\K\= dt\\ Ps (t) - psAW + dta(t)K(t) 



with ||x(t)|| =x(t) -$*(t), 



(43) 



where T denotes the cost functional and K the control 
field. In this notation, the density matrix has been cast 
into vector form for convenience. Here, we demand that 
the trajectory of the density matrix follows a given path 
in time [ps t d{^)\- The function a(t) is used to avoid an 
unbounded intensity of the control field during the opti- 
mization procedure. 

The problem of optimization can now be stated for- 
mally as follows: 



K* 



argmin (J- [K] ) 

K 



(44) 



P22 



1 

0.9 
0.7 
0.5 



10 15 20 s 




10 15 20 Sl 



(c) 

Re[p 23 (t)] 



(d) 

Im[p 23 (t)] 




10 15 20 s 



10 15 20 sc 

FIG. 2: Red (thick solid) lines indicate the desired values 
of the corresponding density matrix elements, (a) The tra- 
jectory of the diagonal element P22, (b) the diagonal den- 
sity matrix element P33 = 1 — P22, (c) the real part of the 
off-diagonal density matrix element Re [P23] = Re[p32], and 

(d) the corresponding imaginary part Im [P23] = — Im[/>32]. 

(e) The concurrence vs time and (f) the control field J(t) vs 
time. 



where K* is the solution of the minimization problem. In 



general, it is quite difficult to find a solution to Eq. (44 1 
because the cost functional can be very complicated as is 
indicated by the structure of Eq. (43 ). We used a Polak- 



Ribicre-type conjugate gradient method to search for a 
minimum of the cost functional. - 

To apply the CG algorithm, the gradient of the cost 
functional has to be calculated. That can be done, e.g., 
using a Hamiltonian approach which leads to a set of con- 
jugate differential equations to those for ps(i)P^ These 
have to be solved backward in time to get the desired 
gradient. A simpler and faster (but numerically less sta- 
ble) method is to parametrize the control field in time 
and to calculate the gradient by a difference method, 

J(t) -» J(t n ) = J„, t n = nh, ne [0,7V] 

N = -f-, ft,... grid spacing 
h 

5T_ _ F{Jg, J„ + AJ„, ...) - J" (Jo, J n , ...) 

SJ n AJ„ 



8T ( 5T 5T 8T 



8 J V(5Jq'(5Ji' ' 8 Jn 



(45) 



which we use here. 

Utilizing the CG algorithm it is possible to optimize 
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FIG. 3: Red (thick solid) lines indicate the desired values of 
the corresponding density matrix elements. (+, blue) data 
points denote the optimized case and (o, black) ones denote 
the initial data or guess, (a) The trajectory of the off-diagonal 
element Re [P23], (b) the imaginary part of pg2, (c) the diago- 
nal element p^g — 1 — P22, (d) the concurrence vs time, (e) the 
control field J(t), and (f) the purity of the system vs time. 



target (trapping) state into the ground state has been 
tested within a Markovian quantum master equation ap- 
proach for a spin-boson coupling between the double 
qubit and its environment. This analysis has shown 
that dcphasing due to the environment modeled by two 
uncorrelated baths of harmonic oscillators can be sup- 
pressed. When the strategy of making the target state 
to the ground state of the system is applied to an ar- 
bitrary maximally entangled state (Bell state), it has 
been shown that a detailed controllability of the (sym- 
metry) of the interdot coupling is required. However, 
an effective isotropic Heisenberg interaction between the 
two spins allows preparation of the selected Bell state 
|^±) = -L(|l) ® |0) ± |0) O |1)). Numerical results 
for the case of steering and trapping into this entan- 
gled state were given. An analysis of the dependence 
of the optimization technique on the temperature of the 
baths has shown that stable entanglement with respect 
to pure dephasing using experimentally realistic qubit- 
qubit couplings is possible over a broad temperature 
regime (< 15 K). Furthermore, the control field was 
subjected to a CG optimization routine which leads to 
improved results: more rapid driving into the target Bell 
state and improved final purity and concurrence. 
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the parameters J n [see Fig. |3je)]. A rapid initial and 
small final increase of the concurrence as well as a faster 
convergence of the density matrix elements were achieved 
by increasing the amplitude of the interdot coupling J(t) 
at intermediate times [see Figs. |3ja)-[3jd)]. It is interest- 
ing that the initial increase of concurrence was gathered 
at the expense of purity at intermediate times [Fig. [3jf )] . 
However, final purity is significantly improved over the 
result from the initial guess for J(t). 



APPENDIX A: CALCULATION OF THE 
CORRELATION FUNCTIONS 



The bath operators of Eq. ( 28 ) in the interaction pic- 
ture read 

r<(*) = (e-^^CO) + e*4°*6»t (0) ). (Al) 



V 


huj sc 


T t f 


— tf,sc/w S c huj c 




M 2) 


0.08 


lmeV 


50 mK 


2.56 ps 20 meV 


1 T 


1 T 



TABLE III: Physical quantities used to calculate the data 
shown in Figs. [3ja,)-j3jf ) . 



Because p^ [see Eq. (30 1] is diagonal in the bosonic oc- 
cupation numbers, we get 

f%)f = tr R {f (i) (*)r w C0pj?} = (A2) 

- E(^0 2 [ e ^ ,(t " t ' ) (^ ) + 1 ) + e ^ ,(t " t ' ) ^ 



IV. CONCLUSION AND OUTLOOK 

It has been shown that entanglement of two electron 
spins in a semiconductor double quantum dot can be 
established and preserved by means of a tunable inter- 
dot voltage barrier. The strategy of transforming the 



where n k denotes the expectation value of the particle 
nu mber operator b^b^\ Now, define the spectral den- 



sity 



27 28 



S[(j- 



(A3) 
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By using Eq. ( A3 1 , we can reformulate Eq. ( A2 ) as 



(f i (t)f i (*')) B = I /wM[(s (t) H + i)e- uM) + 



+ n^(uj)e luj{t ' t ' ) 



(A4) 



If we use the continuum limit [Eq. (29l] and 

1 



n«(w) 



(3 = 



i 



(A5) 
(A6) 



we can write Eq. (|A4| as 



2 a/ 



1 



-iu)(t-t') 



j(t-t') 



This integration can be done analytically yielding 

Eq. m. 



APPENDIX B: DEPENDENCE OF THE 
DYNAMICS ON THE COUPLING STRENGTH r) 

In Figs. |i|a)-|4jc) the oscillations of the diagonal den- 
sity matrix element /?22 with respect to a variation of the 
coupling strength 77 are shown. The control field J(t) is 
shown in Fig. |4jd). As can be seen, there is an anal- 




10 15 20 



10 15 20 s ' 



FIG. 4: Red (thick solid) lines indicate the desired values 
of the density matrix element. The diagonal density matrix 
element P22 vs time are shown for (a) 77 = 0.02, (b) 0.05, 
and (c) 0.10. The control field J(t) is given in (d). For the 
(A7) rema i mn g physical quantities, we used the same values as in 



Table M 



ogy to classical mechanics. For 77 = 0.02 and rj = 0.05, 
the system is "underdamped" , whereas for 77 = 0.10 the 
system is quite close to the "critically damped" case. 
The equilibrium position corresponds to the desired state 
\ip + )(t/j + \ [see Eq. (42 1]. Unfortunately, we cannot ap- 



proach the "overdamped" parameter regime because of 
positivity problems. These problems, e.g., P(t) > 1, 
result from straining the Born-Markov approximation, 
i.e., the system-environment interaction becomes too 
strong .ESI 
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